Helimagnons in a chiral ground state of the pyrochlore antiferromagnets 
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The Goldstone mode in a helical magnetic phase, also known as the helimagnon, is a propagating 
mode with a highly anisotropic dispersion relation. Here we study theoretically the helimagnon 
excitations in a complex chiral ground state of pyrochlore antiferromagnets such as spinel CdC^CU 
and itinerant magnet YMn2. We show that the effective theory of the soft modes in the helical state 
possesses a symmetry similar to that of smectic liquid crystals. We compute the low-energy spin- 
wave spectrum based on a microscopic spin Hamiltonian and find a dispersion relation characteristic 
of the helimagnons. By performing dynamics simulations with realistic model parameters, we also 
obtain an overall agreement between experiment and the numerical spin-wave spectrum. Our work 
thus also clarifies the mechanisms that relive the magnetic frustration in these compounds. 



Introduction. The rich phenomenology associated with 
helical spin ordering has attracted considerable attention 
recently. Magnetic spirals have been shown to play a cru- 
cial role in inducing the spontaneous electric polarization 
in multifcrroic materials [l|. The soft magnetic fluctua- 
tions in a helical spin-density- wave state of itinerant mag- 
net MnSi also holds the key to its non-Fermi liquid be- 
havior ■ Contrary to most well known long-range spin 
orders such as ferromagnetism and antiferromagnetism, 
spin-wave excitations in a helical magnetic order is de- 
scribed by a theory similar to the elasticity of smectic 
liquid crystals [!, HJ. In particular, despite the broken 
symmetry of the helical order is described by a 0(3) order 
parameter, there only exists a single Goldstone mode sim- 
ilar to the smectic-like phonon mode which emerges on 
scales larger than the helical pitch |(| . This low-energy 
gapless excitations, dubbed helimagnons in Ref. 0], ex- 
hibit a highly anisotropic dispersion relation. 

The recent discovery of a novel chiral spin structure in 
spinel CdCr 2 04 has generated much interest both theo- 
retically and experimentally . A very similar copla- 
nar helical order has also been re por ted in the weak itin- 
erant antiferromagnet YMn2 flEl Il6l| . The magnetic ions 
in these compounds form a three-dimensional network 
of corner-sharing tetrahedra (Fig. Q]), known as the py- 
rochlore lattice. The geometrically frustrated spin inter- 
actions in this lattice gives rise to an extensively degen- 
erate ground state at the classical level [l?}- As a result, 
magnetic ordering at low temperatures is determined by 
the dominant residual perturbations in the system. 

The helical magnetic order in CdCr204 is stabilized 
by a combined effect of magnctoclastic coupling and rel- 
ativistic spin-orbit interactions. Despite a rather high 
Curie- Weiss temperature |©cwl ~ 70 K, the spiral mag- 
netic order sets in only at T/v = 7.8 K, indicating a high 
degree of frustration. The magnetic transition at Xjv is 
accompanied by a structural distortion which lowers the 
crystal symmetry from cubic to tetragonal 0) . Recent in- 
frared absorption measurements indicated a broken lat- 
tice inversion symmetry below Tn 12||, consistent with 



the scenario predicted in Ref. [jj that the magnetic frus- 
tration is relieved by the softening of a q = optical 
phonon with odd parity. A long-range Neel order with 
ordering wavevector q = 27r(0, 0, 1) is stabilized by the 
lattice distortion. Moreover, the lack of inversion sym- 
metry endows the crystal structure with a chirality. The 
collinear state is transformed into a magnetic spiral as the 
chirality is transferred to the spins through spin-orbit in- 
teraction. The shifted ordering vector q = 27r(0, 6, 1) is 
consistent with the observed magnetic Bragg peaks. 

In this paper we undertake a theoretical calculation 
of the spin-wave spectrum in the helical ground state 
of pyrochlore antiferromagnet. In particular, our micro- 
scopic calculation goes beyond the phenomenological ap- 
proaches adopted in most studies of magnetic excitations 
in hclimagnets. Based on a well established microscopic 
model for the helical order, we exactly diagonalize the 
linearized Landau-Lifshitz-Gilbert (LLG) equation with 
a large unit cell; the size of the unit cell determines the 
ordering wavevector of the magnetic spiral in the com- 
mensurate limit. By comparing results obtained from dif- 
ferent helical pitchs, we find a spin-wave spectrum that 
is characteristic of helimagnons. Our work thus provides 
a solid example of helimagnons in frustrated pyrochlore 
lattice. We also performed dynamical LLG simulations 
on large finite lattices and found an overall agreement 
with the experimental data on CdCr204. 

Model Hamiltonian. Our starting point is a classical 
spin Hamiltonian that includes nearest-neighbor (NN) 
exchange interactions and spin-orbit coupling manifested 
as the Dzyaloshinskii-Moriya (DM) interaction [3 El 
between NN spins: 
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Kij) Sj • Sj + Dy ■ (Sj x Sj) 



(1) 



The isotropic NN exchange can be recast into 
(J/2) J2 r |M(r)| up to a constant, where M(r) denotes 
the vector sum of spins on a tetrahedron at r. Minimiza- 
tion of this term requires the vanishing of total mag- 
netic moment on every tetrahedron but leaves an exten- 
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sively degenerate manifold. The lattice distortion intro- 
duces the exchange anisotropy K+j through magnetoe- 
lastic coupling. Here we consider tetragonal lattice dis- 
tortions that preserve the lattice translational symmetry. 
In particular, the odd-parity distortion breaks the inver- 
sion symmetry and stabilizes a collinear spin order with 
wavevector q = 2tt(0, 0, 1) HO, ED]. The DM term is al- 
lowed on the ideal pyrochlore lattice, where the bonds are 
not centrosymmetric as required by the so-called Moriya 
rules . Moreover, the high symmetry of the pyrochlore 
lattice completely determines the orientations of the DM 
vectors up to a multiplicative factor. The explicit form 
of the vectors D; 



22, 2 



can be found in Ref. 
To understand how the coplanar helical order is stabi- 
lized, we briefly review the continuum theory based on a 
gradient expansion of the order parameters [7| . The anti- 
ferromagnetic order on the pyrochlore lattice is character- 
ized by three staggered order parameters Lj (i = 1, 2, 3). 
For example, Li = (So + Si — S2 — Sa)/4 describes the 
staggered spins on two opposite [Oil] and [Oil] bonds, 
here the subscript ^ = 0^3 denotes the four sublattices 
of the pyrochlore lattice (Fig.[T|). The order parameter for 
the collinear order is Li(r) = rii e lq r , and L2 = L3 ~ 0, 
where q = 27r(0, 0, 1), and rii is an arbitrary unit vector. 

Assuming a slow variation of staggered magnetizations 
in the spiral, we parameterize the order parameters as: 
L,(r) = e lq ' r (/>i(r) n.;(r), where are three orthogonal 
unit vectors, (pi ps 1 - + </>§), and 4>2,4>3 "C 4>i 0- 
In the J — > 00 limit, vanishing of total magnetization on 
tetrahedra requires 113 || d y hi, and the </> fields can be 
analytically expressed in terms of the director fields n^. 
The effective energy functional in the gradient approxi- 
mation can be solely expressed in terms of the director 
fields rii and 112 || d y iii x fii Q: 

T = D J ' dr hi- (a x d x hi + b x d y hi — c x d z rii) (2) 



-K u \ dr 



(d x ni) 2 + (d v hi) 2 + 2 (<9 2 ni) 2 -(n 2 ■ d z hi) 2 



Here K u and D denote the strengths of exchange 
anisotropy and DM interactions, respectively. The form 
of the DM terms suggests spiral solutions with rii ro- 
tating about one of the principal axes and staying in 
the plane perpendicular to it. These special solutions 
were first obtained in Ref. @- For example, the solu- 
tion rii = (cos Qy, 0, s'mQy) describes a helical order 
with coplanar spins lying in the ac plane; see Fig. [T] 
This spiral magnetic order produces a Bragg scattering 
at Q = 27r(0, 6, 1), consistent with the experimental char- 
acterization of CdCr 2 04 Q. 

The energy functional T actually possesses a contin- 
uous symmetry related to the 0(2) rotational invari- 
ance of the helical axis. To describe the general copla- 
nar spiral solution, we first introduce two orthogonal 
unit vectors e^, which lie in the xy plane. A copla- 
nar spiral rotating about the direction has the form 




FIG. 1: (Color online) Coplanar helical order with spins lying 
in the ac plane and rotating about the b axis. The magnetic 
order corresponds to the solution rii = (cos 6(y), 0, sin Q(y)), 
with 8(y) — Qy + const with Q = 2ir5. It produces a 
Bragg peak Q = 2n(0, S, 1) in neutron scattering, consistent 
with the experimental measurement on CdCr2 04. The ex- 
change anisotropies are -fL[o,i,±i] = K'^ ±1 1] — —2K U , and 

#"[±1,0,1] = K '[o,i,±i] = A '[i.±i,o] = K '[i,±i,o] = K u, where K 
and K 1 denote exchange anisotropy on NN bonds of the green 
and red tetrahedra, and the subscript indicates the bond di- 
rection. 



fii(r) = cos #(£)«;£ + sin 6*(£) c, where £ = r ■ eg . Sub- 
stituting rii(r) into Eq. @ yields an energy density: 
E[6] = -Dd^e + K u {d^ef /4. Minimization of E gives a 
helical wavenumber 



Q = 2nS = d0/d£ = 2D/K U . 



(3) 



This rotational symmetry is accidental and is lifted by 
further neighbor interactions and cubic symmetry field of 
the lattice. A similar 0(3) symmetry also occurs in the 
spiral order of MnSi Q. In addition, all spirals whose 
axis lies in the ab plane are also accidentally degenerate 
with the spiral in which spins rotate about the c axis [3] . 
The experimentally observed (0, S, 1) spiral is selected by 
a large third-neighbor exchange J3 which favors coplanar 
spirals rotating about the a or b axes. The explicit on-site 
magnetizations of a spiral with the pitch vector parallel 
to the 6-axis are 

S M (r) = a sin(Q • r + ip a ) + c cos(Q • r + ip u ), (4) 

where (j, = 0, 1, 2, 3 is the sublattice index, tpo = ip, ipi = 
ir5 + ip, ip2 = n + ip, and </?3 = 7r(l + S) + with ip is an 
arbitrary constant related to the U(l) symmetry of the 
coplanar spiral, i.e. 8(r) = Q • r + ip. 

Helimagnons in the continuum limit. We next consider 
a phcnomcnological approach to low-energy magnetic ex- 
citations in the helical ground state. For a spiral with a 
fixed axis the obvious soft modes are the U(l) phase fluc- 
tuations ip associated with the coplanar director field. A 



3 




4 
3.5 

3 


-2 

xlO 


( r ) 


1 1 
6=0.2 + 
" 5=0.1 
6=0.05 * 


i 

- 


2.5 






Ei 2 






1.5 






1 






0.5 









< i i 


i 



xlO 



0.005 0.01 0.015 
(0 k) 

(f) 




0.04 



0.04 0.08 



0.12 



0.16 



FIG. 2: (Color online) Spin-wave dispersions along three principal axes [panels (a)-(c)] obtained from numerical exact diago- 
nalization of a large unit cell. The results are fitted to the phenomenological dispersion u = \J ak 2 + [3k 4 . The magnon energy 
is measured in unit of JS. Panels (d)-(f) show the dispersion parameters a and /3 as a function of the helical wavenumber S. 
These data are fitted to a =const, a oc 8 2 , and /3 cx S~ 2 , respectively. 



naive guess of the effective energy functional of the soft 
mode would be J 7 cx K u J dr\ V-i/>| 2 , similar to that of pla- 
nar magnets. However, a phase variation ip = a£ corre- 
sponding to a rotation of the helical axis is a zero mode, a 
manifestation of the U(l) symmetry. Yet, it costs a finite 
amount of energy oc K u a 2 . Consequently, there cannot 
be any (d^ijj) 2 term in the effective energy density. The 
correct energy functional for the soft modes is [3|45| 



T=-\dv 



Cz (^) 2 +C||(^) 2 +^(^V) 2 - 



(5) 



where c\\, c± ~ K v 



J3 are elastic constants, and 



the true soft mode is a generalized phase accompanied 
by a rotation of the director field. We have also included 
a zero wavevector ferromagnetic mode m which is soft 
due to spin conservation. It is worth noting that the first 
term in Eq. ([5]) is absent in conventional helical order 
discussed in Ref. 0. As discussed above, it originates 
from a rather large J3 that gives penalty to variations 
along the z direction. Schematically the coarse-grained 
sublatticc magnetizations are 

S M (r) ~ S„(r) + m(r)b (6) 
+ ^(r) [cos(Q • r + ^) a - sin(Q • r + ip^) c] , 

The spin dynamics is governed by a generalized Landau- 



Lifshitz equation [J 
dt 



■7S„(r) x 



SJ 7 



(7) 



where 7 is a constant. Substituting S M and T into the 
Eq. (JT)) yields coupled differential equations dtip = ^rm 
and d t m = — 7(— C||9| — c z d 2 z + c±d^/Q 2 )^. Assuming 
a space-time dependence exp(ik ■ r — iuit), we obtain a 
highly anisotropic dispersion relation [f| 



w(k) = ir 1 ' 2 ^c z k 2 + Cll k 2 +c ± ki/Q 



(8) 



For wavevectors parallel to the spiral and the c axes the 
dispersion is linear as in an antiferromagnet, while it is 
quadratic for wavevectors parallel to the a axis. 

Exact diagonalization. At the microscopic level, the 
crystal symmetry anisotropy is expected to modify the 
helimagnon dispersion ([5]). In addition, the particular 
spiral (0,(5, 1) is stabilized by a large J3 which would in- 
crease the spin stiffness along the c axis • To take into 
account these effects, here we perform exact diagonaliza- 
tions with a large unit cell to investigate the low-energy 
acoustic-like magnons and compare the results with pre- 
dictions based on the helimagnon theory. 

To begin with, we first choose a commensurate unit cell 
of a size which contains 2 x A tetrahedra, corresponding 



4 



to an ordering wavevector q = 27r(0, S, 1) with 6 = 1/A. 
Here the factor of 2 accounts for the staggering of spins 
along the c axis. The DM constant is then determined 
according to the expression D = (3/>/2)K u ta,nn5 [the 
discrete version of Eq. ([3])] for a given fixed K u . Next 
we introduce small perturbations to spins in the ground 
state: Si = S i>i + &i, where i>i is a unit vector par- 
allel to the magnetic moment at site in the ground 
state, and <x; _L i>i denotes transverse deviations. Note 
that because of this constraint, there are only two de- 
grees of freedom associated with cr^ at each site. As- 
suming a time-dependence er^ ~ exp(ik ■ — icot) for 
the eigenmodes, the linearized LLG equation can be cast 
into an eigenvalue problem T(k) <r = —ico^cr, where <r = 
(<7i, (72, • • • , <J n ) is a vector of dimension n = 2x 2A x 4, 
and T is an n x n matrix. We then numerically diago- 
nalize the matrix T(k) and obtain the lowest eigenvalues 
in the vicinity of k = 0, corresponding to the soft-mode 
around the ordering vector of the magnetic spiral. 

The magnon dispersions along the three principal di- 
rections in the vicinity of the ordering wavevector Q = 
27r(0, S, 1) are shown in Fig. HJa)-(c). The number of 
spins in the unit cell corresponding to used spiral pitch 
S is N s = 40, 80, 160. To incorporate the anisotropy 
due to the cubic symmetry of the lattice, the obtained 
spectra arc fitted to the phenomenological dispersion: 
u)(k) = y ak 2 + /3fc 4 . The dispersion along the b and c 
axes can be well approximated by a linear relation, while 
the spectrum along the a-axis exhibits a predominant 
quadratic behavior at small S similar to that of a fer- 
romagnet. The dispersion parameters a and /3 obtained 
from the fitting arc shown in Fig.[^d)-(f) as a function of 
the helical wavenumbcr 5. The large spin-wave velocity 
along the c axis can be attributed to a large J3, consis- 
tent with the analysis of gradient expansion [7|. More 
importantly, the (5-dependence of parameters for disper- 
sion along the a-axis exhibits a scaling relation a ~ S 2 
and P ~ 1/5 2 , reminiscent of that of a helimagnon 

Magnons in CdCr^O^. Although the exact diagonal- 
ization method discussed above provides a rather accu- 
rate magnon spectrum, the complicated band-folding due 
to large unit cell makes it difficult to compare the results 
with experiment. Instead, here we employ a less accu- 
rate but more direct approach based on simula ting the 
real-time LLG equations in a large finite lattice 241 ] : 



dt 



= S, : 



on 



dt 



(9) 



where ac is the dimensionless Gilbert-damping coeffi- 
cient. The dynamics simulation is initiated by a short 
pulse localized at = 0. Since the system is driven by a 
white source with fiat spectrum, we expect magnetic ex- 
citations of various energy and momentum are generated 
in our simulations. The spin-wave spectrum is then ob- 
tained via the Fourier transform of the simulation data. 
The magnon spectrum of CdCr204 has been measured 




(1k0) 

FIG. 3: (Color online) Simulated spectra of (a) Neel state 
with q = 27r(0, 0, 1) and (b) helical order with q shifted to 
27r(0, 5, 1) where the helical pitch S — 0.1 due to a finite D = 
0.14 meV. Other parameters used in the simulation are J = 
1.35 meV, K u = 0.21 meV, and J 3 = 0.28 meV. (b) Fitted 
spin wave spectrum to the neutron scattering measurement 
(white circles) [!|. The finite DM interaction splits each band 
into three as shown in panel (b). The vertical peaks at the 
ordering wavevectors are numerical artifacts. 



using inelastic neutron scattering in Ref. |8J. Here we 
present our numerical calculations and compare the re- 
sults with the experimental findings. The measured heli- 
cal pitch (5 ~ 0.1 can be used to fix the ratio of the DM in- 
teraction to the exchange anisotropy K u . Then the abso- 
lute value of K u can be estimated by the magnon energies 
at the zone center. As for further-neighbor exchanges, ab 
initio calculations find a negligible J2 and a quite large 
J3 with antiferromagnetic sign 0, [25| . Using the follow- 
ing set of parameters: J = 1.35 meV, K u = 0.21 meV, 
D = 0.14 meV, and J3 = 0.28 meV, we obtained very 
good agreement with the experimentally measured spec- 
trum, as shown in Fig- EJb). 

Concluding discussions. We have studied the spin- 
wave excitations in a complex helical magnetic order on 
the pyrochlore lattice. Our exact diagonalization calcu- 
lation of the low-energy dispersion in the commensurate 
limit of the long-period spiral clearly demonstrates the 
existence of helimagnon excitations in the frustrated py- 
rochlore antiferromagnet. This novel magnetic excitation 
can be observed in spinel CdCi'204 and itinerant magnet 
YMn2. We have performed dynamical LLG simulations 
with parameters inferred from structural data and ab ini- 
tio calculations. The numerical magnon spectrum agrees 
well with the measured dispersion in CdCr204. Our work 
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thus helps clarify the underlying mechanisms that stabi- 
lize the observed helical order and provides a quantitative 
description of the model Hamiltonian. 
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